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ABSTRACT 

We perform a multi-wavelength polarimetric study of the quasar GTA 102 during an extraordinarily 
bright 7 -ray outburst detected by the Fermi Large Area Telescope in September-October 2012 when 
the source reached a flux of F>ioo MeV = 5.2 ± 0.4 x 10“® photons cm“^ s“®. At the same time the 
source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec- 
scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity 
Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 
2014. We And that the y-ray outburst is coincident with flares at all the other frequencies and 
is related to the passage of a new superluminal knot through the radio core. The powerful y-ray 
emission is associated with a change in direction of the jet, which became oriented more closely to 
our line of sight {6 ~1.2°) during the ejection of the knot and the y-ray outburst. During the flare, 
the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, 
which we associate with the path followed by the knot as it moves along helical magnetic field lines, 
although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We 
locate the y-ray outburst a short distance downstream of the radio core, parsecs from the black hole. 

This suggests that synchrotron self-Gompton scattering of near-infrared to ultraviolet photons is the 
probable mechanism for the y-ray production. 

Subject headings: galaxies: active — galaxies: jet — (galaxies:) quasars: individual (CTA102) — 
techniques: interferometric — techniques: photometric — techniques: polarimetric 


1. INTRODUCTION 

The Blazar CTA 102 (B2230-I-114) is classified as a 
highly polarized quasar (HP Q), characterized by optica l 


polari 
Its hi 

zation exceeding 3% (|Moore & Stockman 19811). 

gh variabilitv at optical wavelengths (fPica et al. 

1988 

Osterman Meyer et al.|2009 

) and its spectral prop- 


erties identity it also as an optically violent variable 
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(OVV) quasar ( 

Maraschi et al. 

1986 

). Microvariability 

of GTA 102 at o 

ptical waveiengt 

hs has been investigated 

by Osterman Meyer et al. ( 

2009 

), who found that faster 


The variability in this source occurs at other frequen¬ 
cies as well: flares at cm and mm wavelengths have been 
registered in the past, as well as an X-ray flare d etected 


by the Rossi X-ray Tim ing Explorer in late 2005 (Oster- 


man Meyer et al.|2009|. A radio flare in 1997 was related 


to the ejection or a new knot from the core (ISavolainen 


et al. 2002 Rantakyro et al. 2003 Jorstad et ai. 2005 f 


and another, in 2006, was connected with both the ejec¬ 
tion of a new superluminal feature and the interaction 
betw een this component and a recollimation shock at 0.1 
mas (|Fromm et al.||2011|). 

The radio morphology on arcsecond scales (from im¬ 
ages with the Very Large Array at 6 and 2 cm) reveals a 
central co re with two weak components located at oppo - 


iStanghellini et al.||1998|). 
TA 102 has been regularly 


site sides (Spencer et al. 1989 

At higher angular resolution, O'iA 102 has been regularly 
observed since 199 5 within the VLBA 2cm-Survey (e.g. 
Zensus et al.||2002') and its successor, the MOJAVE pro- 


gram (e.g., Lister et al. 2009). MOJAVE images show 
that the jet in O'iA 102 extends toward the southeast 
up to ^15 mas from the core, which corresponds to a 
de-projected distance of ^ 2.7 k pc using the estimate d 
viewing angle of 2.6° obtained by|Jorstad et al.|(|2005D. 
Kinematic studies of the MOJAVE dat a report ap- 


parent velocities between 1.39c and 8.64c (Lister et al. 
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20131. Higher apparent speeds, up to Papp ~18c, have 
been reported in higher-resolution Very Long Baseline 


Array (VLB A) observations at 43 GHz by Jorstad et al. 
(120011120051) . Apart from superluminal features, the jet 


ot Cl A 102 also displays standing features: two sationary 
components, Al and C, have been observed at a distance 
of ^ 0.1 and 2 m as from the core, respectively (Jorstad 


et al.|200l 2005|), and interpreted as recollimation shocks 
(I'romm et al.||2013b I. 

Recen t MOJAVE polarime tric multifrequency obser¬ 
vations (Hovatta et al. 2012) reveal a rotation measure 
gradient across the jet width at 7 mas from the core, 
which suggests the presence of a helical magnetic field 


in the j et. A similar result is reported in Fromm et al. 
(2013b), where the different observing frequencies reveal 


bends and helical structures in many parts of the jet. 

CTA 102 was detected by the Fermi Gamma-ray Space 
Telescope in the hrst Large Area Telescope (LAT) catalog 
with a flux (E > 100 MeV) of 14. 70 ± 0.97 x 10"® pho¬ 
tons cm“^ s“^ ( Abdo et al.| 2010 '^, and confirmed in the 
second catalog ( Ackermann et al.|2011 ). In late 2012 the 
blazar exhibited an extraordinarily bright y-ray outburst 
detected by the LAT, reaching a flux of 5.17±0.44x 10“® 
photons cm“^ between 0.1 and 200 GeV. During 
the 2012 event, an unprecedented optical and NIR out¬ 
bur st was observed, as rep orted by Larionov et al. (2012 1 
and Garrasco et al. (2012), respectively. 

In this paper we perform a radio to y-ray multi¬ 
wavelength analysis in order to study the flaring activity 
of CTA 102 during the 2012 event. In Section 2, we 
present the multi-wavelength data set collected for the 
analysis, and we describe the methods used to reduce 
the data. In Section 3 we describe the characteristics of 
the emission at the different frequencies during the 2012 
flare event. In Section 4 we study the kinematics and the 
flux density variability of the parsec scale jet. In Section 
5 we perform the discrete cross-correlation analysis be¬ 
tween light-curves. In Section 6 we analyze the polarized 
emission of the source at mm and optical wavebands. We 
present our discussions and conclusions in Sections 7 and 


(2) 230 GHz (1.3 mm) and 86.24 GHz (3.5 mm) with 
the 30 m Telescope of Institut de Radioastronomie 
Millimetrique (IRAM) at the Pico Veleta Observatory 
(Spain); (3) 43 GHz (7 mm) with the Very Long Base¬ 
line Array (VLBA); and (5) 37 GHz (8 mm) with the 13.7 
m Telescope at Metsahovi Radio Observatory of Aalto 
University (Finland). 

Near-infrared photometric data (JHK filters) were ob¬ 
tained at the Perkins Telescope at Lowell Ob servatory 
(Flagstaff, AZ) using the MIMIR instrument (Glemens 


et al. 2007| and at the 1.1m Telescope of the Mam Astro¬ 
nomical Observatory of the Russian Acade my of Sciences 


located at Campo Imperatore, Italy (see Hagen-Thorn 
et al. [20(38 for details). 


We have collected optical photometric data in UBVRI 
bands from numerous telescopes: (1) the 2.2 m Telescope 
of Calar Alto Observatory (Almen'a, Spainj^ (2) the 
2 m Liverpool Telescope of the Observatorio ael Roque 
de Los Muchachos (Canary Island, Spain); (3) the 1.83 m 
Perkins Telescope of Lowell Observatory (Flagstaff, AZ); 

(4) the 1.54 m and 2.3 m telescopes of Steward Ob¬ 
servatory (Mt. Bigelow and Kitt Peak, AZj^ (5) the 
40-cm LX-200 Telescope of St. Petersburg slate Uni¬ 
versity (St. Petersburg, Russia); ( 6 ) the 70 cm AZT- 
8 Telescope of the Crimean Astrophysical Observatory 
(Nauchnij, Ukraine); (7) the 1.5 m Kanata Telescope in 
Higashi-Hiroshima Observatory (Japanl^ and ( 8 ) the 
Ultraviolet and Optical Telescope (UVO^ on board the 
Swift satellite. Optical data are in part also in polari¬ 
metric mode, mostly in R band, except for items (4) and 

(5) listed a bove (see Schmidt et al. 1992 Hagen-Thorn 
et al.|2008| respectively, tor details). In the U V range we 
use UVO'i' data from Swift in the three available filters: 
UVW 2 (2030 A), UVM 2 (2231 A) and UVWl (2634 A). 
We have also obtained X-ray data in the energy range 
0.3-10 keV from the X-Ray Telescope (XRT) on board 
of Swift satellite. 

At the highest photon energies considered here, we 
have analyzed y-ray data from the Large Area Telescope 
(LAT) of the Fermi Gamma-ray Space Telescope. 


We adopt the cosmologic al values from the most re¬ 


cent Planck satellite results (Planck Collaboration et al. 
20141: 0.3, Ua= 0.7, and FIq = 68 km s“^ Mpc“^ 


With these values, at the redshift of CTA 102 ( 2 ;=!.037), 
1 mas corresponds to a linear distance of 8.31 pc, and a 
proper motion of 1 mas yr“^ corresponds to an apparent 
speed of 55.2c. 

2. OBSERVATIONS AND DATA REDUCTION 


We are interested in studying the multi-spectral be¬ 
havior of CTA 102 during the 7 -ray flare that occurred 
between 2012 September and October. For this, we have 
collected data from millimeter to 7 -ray wavelengths, ex¬ 
tending our study from 2004 June to 2014 June. In par¬ 
ticular, the Fermi data extend from the start of the mis¬ 
sion (2008 August) to 2013 September, X-ray and UV 
data cover the observing period from 2009 August to 
2013 June, optical and NIR data from 2004 June to 2013 
October, and the combined radio light curves cover the 
entire period from 2004 June to 2014 June. 

In the mm-wave range, we use data at (1) 350 GHz 
(0.85 mm) and 230 GHz (1.3 mm), obtained with the 
Submillimeter Array (SMA) at Mauna Kea, Hawaii; 


2 . 1 . 'y-ray Data Analysis 

We have analyzed the 7 -ray data of the field contain¬ 
ing CTA 102 from the Fermi LAT from 2008 August 
to 2013 September, producing a light curve between 0.1 
and 200 GeV with an integration time of 1 day. We 
used the Fermi Science Tools version v9r33p0 and in¬ 
strument response function P7REP_SOURCE_V15, con¬ 
sidering data inside a region of interest (ROI) of 15° 
radius centered on CTA 102. We followed the un¬ 
binned likelihood procedure as described in web page 
http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/. 

The procedure starts with the selection of good data 
and time intervals through the tasks gtselect and gtmk- 
time, and follows with the creation of an exposure map 
for each day (tasks gtltcube, gtexpmap) and the modeling 
of data through a maximum-likelihood method (gtlike). 
In this last step, we used a model that includes CTA 102 


Observations performed under the MAPCAT (Monito ring 
AGN with the Calar Alto Telescopes), see |Agudo et al! 

Data taken from t he Steward Observatory monitoring project, 
see [Smith et al~ ]l |200^ 


Data published in 


Itoh et al. 


120131 
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Fig. 1.— Fermi LAT 7 -ray light curve between 0.1 and 200 GeV, 
with an integration time of 1 day. Red points represent the detec¬ 
tions (TS>10) and gray arrows correspond to upper limits when 
the source is not detected (TS<10). 

and 42 other point sources inside the ROI, accor ding to 


the second Fermi Gamma-ray Catalog, (2FGL; Acker- 
mann et al. 20111. The model also takes into account 
the dittuse emission from our Galaxy {glLiem-v05.fit) and 
the extragalactic 7 -ray background {iso_source-v05.txt). 
We searched for the flux normalization of CTA 102 by 
fixing the spectral index of the other point sources while 
varying both the flux and spectral index of our target. 
The 7 -ray spectrum of CTA 102 was modeled with a log 
parabola curve corresponding to the spectral model given 
in the 2FGL catalog. We considered a successful detec¬ 
tion when the test statistic TS > 10 , which corres ponds 
to a signal-to-noise ratio > 3-cr (Nolan et al.||2012 1 . 


2 . 2 . X-ray, UV, Optical and Near Infrared 

We collected X-ray and UV data from 2009 August to 
2013 June from the Swift archive. The X-ray data in the 
energy range 0.3-10 ke V were calibrated followi ng the 

([Mil). 

U 


procedure described in Williamson et al. 
UVOT data reduction was performed via the 


The 

VOT- 


SOURCE tool, with a sigma value of 5 adopted to com¬ 
pute the background limit. Optical and NIR data were 
reduced and calibrated following the procedures outlined 


Jorstad et al. (2010). 

All of the magnitudes of the optical and NIR data have 
been corrected for Galactic extinction with values re¬ 
ported in the NASA Extragalactic Da tabase (NEPf^ 
for each filter (Schlafly & Finkbeiner 2011). For tfe* 
UV data, we obtained the absolute extinction value at 


each w avelength A(A) from equation (1) in Gardelli et al. 
(1989). After the correction we transforme d magnitudes 
into fluxes using the formula reported in Mead et al. 


(1990) and Poole et al. (2008). 


2.3. Photo-Polarimetric Millimeter VLBA and Single 
Dish Data 

http://ned.ipac.caltech.edu/ 


Multi-epoch very long baseline interferometer (VLBI) 
images provides ultra-high angular resolution that can 
be used to determine the location in the jet where flar¬ 
ing activity occurs. We therefore have collected data 
from the VLBA-BU-BLAZAR prograirp^ which consists 
of monthly monitoring of 7 -ray bright olazars with the 
VLBA at 43 GHz (7 mm). The dataset consists of 80 
total and polarized intensity images from 2007 June to 
2014 June. We restore the images with a common con¬ 
volving beam of 0.4x0.2 milliarcseconds (mas). Since 
the resolution of the longest baselines of the VLBA is 
less than half of these dimensions, we employ model fit¬ 
ting to define and determine the parameters of the very 
fine-scale structure. The data reduction was performed 
with a combination of the Astronomical Image Process- 


(Difmap), as described in 

Jorstad et al. 

(2005 

). The 

electric vector position angle (EVFA) calibration follows 

the procedure discussed in 

Jorstad et al. 

2005 

, which 


combines the comparison between VLA and VLBA inte¬ 
grated EVPA values at those ep ochs for which VLA data 
are available with the method of Gomez et al. (2002) that 
utilizes the stability of the instrumental polarization (D- 
terms). 

The IRAM 30m Telescope’s total flux and polarimet- 
ric data in this paper were acquired under the POLAMI 
(Polarimetric AGN Monitoring at the IRAM 30m Tele¬ 
scope) program (see Agudo et al. MNRAS in prep.) and 
reduced and calibrated following the pr ocedures intro¬ 
duced in|Agudo et al.l(|2006l|2010l |2014 ). 

The Submillimiter Array data of G'i'A 102 came from 
an ongoing monitoring program at the SMA to deter¬ 
mine the fluxes of compact extragalactic radio so urces 
that can be used as calibrators at mm wavelengths (Gur- 
well et al. ]2007t . Observations of available potential cal¬ 
ibrators are from time to time observed for 3 to 5 min¬ 
utes, and the measured source signal strength calibrated 
against known standards, typically solar system objects 
(Titan, Uranus, Neptune, or Gallisto). Data from this 
program are updated regularly and are available at the 
SMA websitePH 

3. MULTI-WAVELENGTH OUTBURST 

Figure [T] displays the 7 -ray light curve of GTA 102 
in the energy range 0.1-200 GeV during the period of 
major activity (2011 June-2013 Ap ril) obtained with an 
integration time of 1 day. Following Jorstad et al. (2013) 
we can define a 7 -ray outburst as a period when the 
flux exceeds a threshold of 2 x 10 “® photons cm“^ s“^. 
Although this is an arbitrary limit, it conforms to a visual 
inspection of the 7 -ray light curve of GTA 102. 

The first outburst takes place in 2011 June (MJD 
55719-55721), when the source displays a one-day peak 
flux of 3.1 ± 0.37 X 10“® photons cm“^ s“^. The second, 
brightest outburst occurs at the end of 2012 September 
(2012.73), when the source remains above 2 x 10“® pho¬ 
tons cm“^ s“^ for 14 days (from MJD 56188 to 56202), 
reaching a peak of 5.2 ± 0.4 x 10“® photons cm“^ s“^ on 
MJD 56193. During this outburst the 7 -ray flux increases 
by a factor of ten in just 6 days. The third flare occurs 
in April 2013 (MJD 56387-56394) and lasts 8 days. On 

http://www.bu.edu/blazars/research.html 
http://smal.sma.hawaii.edu/callist/callist.html 
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Fig. 2.— Light curves of CTA 102 from 7 -ray to millimeter wavelengths. From top to bottom: 7 -ray, X-ray, UV, optical, NIR and 
millimeter-wave data. Left panel: Data from 2004 May to 2014 January. Right panel: Expanded view during the 7 -ray outburst between 
2012 August and November. 


this occasion (MJD 56392), the blazar reaches a peak of 
2.9 ± 0.4 X 10“® photons cni“^ s“^. 

We compare the daily 7 -ray light-curve with the X-ray, 
UV, optical, NIR, and radio light-curves in Figurej^ Ta- 
ble[l]lists the multi-wavelength data used in our armlysis. 
Thebrightest 7 -ray outburst, in 2012 September, is ac¬ 
companied by similarly bright flares at all of the other 
wave bands. This is, however, not the case for the other 
two orphan 7 -ray flares ( 7 -ray outbursts with no corre¬ 
spondence at any of the other observing bands), with the 
exception of a mm-wave flare that follows the third 7 -ray 
flare in 2013 April. 

Analyzing the multi-wavelength flare in 2012, we ob¬ 
serve that the X-ray light curve contains a double-peak 
structure, where the first peak is almost coincident with 
the 7 -ray outburst and the second peak occurs ~50 days 
later. The limited sampling of the X-rays prevents a 
deeper analysis of the overall X-ray behavior associated 
with this flare. The UV and optical bands exhibit a 


rapid and pronounced increase in the light curves coin¬ 
ciding with the 7 -ray flare. In the optical light curve, 
we distinguish a secondary, weaker flare after ~50 days, 
close to the second X-ray peak, as well as a third, smaller 
outburst that occurs ^70 days after the second peak. In 
the UV it is also possible to distinguish a secondary flare 
delayed by ^50 days with respect to the main flare, but 
the sampling of the data is insufficient to specify the be¬ 
havior in more detail. In the NIR light curve, we observe 
a large flare coincident with the 7 -ray flare, but there is 
no further sampling after this. A detailed analysis of the 
NIR flare shown in Figure reveals that the event con¬ 
sists of three sub-flares covering almost the entire period 
of high 7 -ray flux from MJD 56193 to 56202. 

The radio light curve also exhibits an increase in flux 
density during the 2012 7 -ray outburst, but with a much 
longer time scale, lasting ~200 days. The 1 mm light 
curve and 7 mm light curves of features CO and Cl peak 
on ^ MJD 56230, about 1 month after the 7 -ray flare. 
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TABLE 1 

Multi wavelength Data 

7 -ray data 



Epoch 

Flux 

Energy band 

(MJD) 

(phot / cm^ / sec) 

(GeV) 

54684.2 

2.15e-07±l.lle-07 

0 . 1-200 

54688.2 

5.68e-07±2.30e-07 

0 . 1-200 

54689.2 

10.00e-08±8.30e-08 

0 . 1-200 


X-ray data 

Epoch 

Flux 

Energy band 

(MJD) 

(erg/cm^/sec) 

(KeV) 

53509.4 

5.14e-12±4.74e-13 

0.3-10 

54210.2 

4.07e-12±4.13e-13 

0.3-10 

54212.1 

5.97e-12±7.67e-13 

0.3-10 


(This table is available in its entirety in machine-readable form.) 



Fig. 3.— An expanded view of the NIR light-curve in the J, H, 
and K bands during the period of the main 7 -ray outburst. 

The 3 mm light curve follows a similar trend, starting to 
increase on ~ MJD 56000. Our limited time sampling 
between MJD 56208 and 56412 shows a peak on MJD 
56207, very close to the 7 -ray flare, although we cannot 
rule out the possibility that the actual peak is closer to 
that at 1 mm. 


4. THE PARSEC-SCALE JET 
4.1. Physical Parameters of Components 

VLBA images of CTA 102 at some selected epochs are 
displayed in Figures [4][^ To carry out an analysis of the 
jet kinematics and nux density variability, we have fit 
with Difmap the complex visibilities with a model source 
consisting of components described by circular Gaussian 
brightness distributions. For each epoch, we obtained a 
model-fit that provides information about the flux den¬ 
sity (S), distance (r) and position angle (0) relative to 
the core, and FWHM size (a) of each component. The 
core (labeled CO), considered stationary over the entire 
period, is identified with the unresolved component in the 
north-western (upstream) end of the jet. It is the bright¬ 
est feature in the jet at most of the epochs. Polarization 
information has been obtained with an IDL program that 
calculates the mean values of the degree of polarization 
(m) and electric vector position angle (EVPA, y) over 


the image area defined by the FWHM size of each com¬ 
ponent. The uncertainties of both m and y correspond to 
the standard deviations of their respective distributions. 
Model-fit parameters for all components and epochs are 
reported in Table 

The accuracy of the model-fit parameters for each 
component depends on its brightness temperature, so 
that smaller uncertainties are expected for more compact 
components and higher flux densities. We have therefore 
established a criterion for quantifying the errors in the 
model ht parameters that is directly related to the ob- 
served brightness t emperature, Th = 7.5 x 10®S'/a^ (e.g., 
Jorstad et al. 2005), where Tf, is measured in Kelvins (K), 
i>' in Janskys (Jy), and a in milliarcseconds (mas). 

First, we select a representative sample of epochs and 
components with a wide range of 5, a and r values. For 
each one of these components, we compute the error in 
the fitted parameters by analyzing how the reduced y^ 
of the fit and resulting residual map change when vary¬ 
ing the fitted parameters one at a time. We set a limit 
on the maximum allowed variation of the reduced y^ of 
20%, corresponding to an increase by a factor of ~1.5 
in the peak levels of the residual map. According to 
this criterion, we assign a series of uncertainties in posi¬ 
tion and flux density to the components in the sample. 
We then relate the derived uncertainties with the mea¬ 
sured brightness temperature, obtaining the following re¬ 
lations: 

« 1.3x lO^T-O-e, (1) 


as ~ 0.09 T. 


- 0.1 


( 2 ) 


where a^y and as are the uncertainties in the position 
(right ascension or declination) and flux density, respec¬ 
tively. These relations have been used to compute the 
errors in the position and flux density for all of the fit¬ 
ted components. To account for the errors in the flux 
calibration, we have added in quadrature a 5% error to 
the uncertainty in flux density. The uncertainties in the 
sizes of components are also expect ed to depend on thei r 
brightness temperatures. Following Jorstad et al. (2005), 
we have assigned a 5% error to the sizes of the rnajority 
of components (those with flux densities > 50 mJy and 
sizes of 0.1-0.3 mas) and a 10% error for more diffuse 
components. 

Plots of separation and flux density versus time for the 
model-ht jet components, including the core, are pre¬ 
sented in Figures and Besides the core, we have 
identified seven main components that could be traced 
reliably across multiple epochs. Component El, located 
at ~2 mas from the core, is a weak and extended feature 
that appears to be quasi-stationary across some epochs, 
or to move with a significantly slower velocity than other 
moving components (see Figure [ 7 ]). A stronger and more 
compact component, Cl, can be distinguished from the 
core at most of the observed epochs, located at a mean 
distance r ~0.1 mas. Both qu asi-stationary features have 
been observed previou sly by Jorstad et al. (20051 and 
Fromm et al. (2013b|a), and interpreted as recolliniation 
shocks m the jet. We .dentify hve other moving compo¬ 
nents, Nl, N2, N3, N4 and SI. Component SI seems, 
however, to have a different nature: it appears to form 
in the wake of component Nl at a distance of ^0.5 mas, 
and it is observed over only four epochs afterwards. Its 
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Relative Right Ascension (mas) 


Linearly Polarized Intensity (mJy/beam) 
3.27 38.63 73.99 109.35 144.72 



Fig. 4.— Sequence of 43 GHz VLBA images displaying epochs from 2009 February to 2010 August, when we observe the appearance of 
component Nl. The images are restored with a common beam of 0.4x0.2 mas at —10° and are separated by a distance proportional to 
the time elapsed between observing epochs. Left panel: contours (total intensity) are traced at 0.003, 0.008, 0.04, 0.1, 0.3, 0.6, 1.2, 1.8, 2.5 
and 3.0 Jy/beam and 7pea,fc=4.2 Jy/beam. Red circles represent modelfit components. Right panel: same contours (total intensity) as in 
left panel plus colors that represent linearly polarized intensity and white sticks symbolizing linear polarization angle. 
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Fig. 5.— Sequence of total intensity 43 GHz VLBA images from 2012 October to 2014 June, covering the epochs from which we start 
observing component N4. Peak intensity is -fpea,fc=3.6 Jy/beam and contours are traced at 0.003, 0.008, 0.04, 0.1, 0.3, 0.6, 1.2, 1.8, 2.5 and 
3.0 Jy/beam. 


Linearly Polarized Intensity (mJy/beam) 

1.47 12.14 22.81 33.48 44.15 



Fig. 6.— Same as Figure[^with linearly polarized intensity in colors and white sticks symbolizing linear polarization angle. 




TABLE 2 

VLBA 43 GHz Model-fit Components’ Parameters 


Epoch 

(year) 

Epoch 

(MJD) 

Elux 

(mjy) 

Distance from 
CO (mas) 

Pos. Angle 

n 

Major 
axis (mas) 

Degree of 
polarization (%) 

EVPAs 

n 

Component CO 

2007.45 

54264.5 

3086 ± 313 

- 

- 

0.017 ± 0.001 

- 

- 

2007.53 

54294.5 

3423 ± 347 

- 

- 

0.034 ± 0.002 

1.5 ± 0.2 

71.2 ± 5.7 

2007.59 

54318.5 

2340 ± 239 

- 

- 

0.038 ± 0.002 

2.4 ± 0.1 

76.4 ± 5.4 

2007.66 

54342.5 

3163 ± 321 

- 

- 

0.045 ± 0.002 

1.6 ± 0.1 

71.9 ± 6.6 

2007.74 

54372.5 

2743 ± 279 

- 

- 

0.045 ± 0.002 

1.4 ± 0.1 

82.9 ± 8.4 


Component Cl 

2007.66 

54342.5 

296 ± 37 

0.07 ± 0.01 

119.2 ± 3.2 

0.062 ± 0.003 

1.6 ± 0.1 

65.5 ± 7.1 

2007.74 

54372.5 

307 ± 37 

0.09 ± 0.01 

104.4 ± 4.8 

0.031 ± 0.002 

2.2 ± 0.1 

-80.8 ± 6.3 

2007.83 

54405.5 

168 ± 25 

0.13 ± 0.02 

119.5 ± 3.2 

0.116 ± 0.006 

2.6 ± 0.2 

79.3 ± 5.9 

2008.04 

54482.5 

575 ± 64 

0.08 ± 0.01 

146.2 ± 1.9 

0.070 ± 0.003 

0.7 ± 0.1 

37.8 ± 5.5 

2008.16 

54524.5 

894 ± 96 

0.09 ± 0.01 

-166.3 ± 7.4 

0.062 ± 0.003 

0.5 ± 0.1 

38.2 ± 5.6 



(This 

table is available 

in its entirety in 

machine-readable form.) 



TABLE 3 

Kinematics of Moving Jet Features 


2007 2008 


2009 
' Ini 


2010 


Name 

N.Epoch 

(mas/yr) 

l^app 

(c) 

T 

J- ej 

(year) 

Nl 

26 

0.27±0.01 

14.9±0.2 

2009.12±0.02 

N2 

18 

0.35±0.01 

19.4±0.8 

2010.65±0.07 

N3 

10 

0.49±0.03 

26.9±1.8 

2011.96±0.07 

N4 

6 

0 . 21 ± 0.02 

11.3±1.2 

2012.49±0.11 


2007 2008 2009 2010 2011 2012 2013 2014 



Fig. 7.— Distance from the core vs. time for the 43 GHz model- 
fit components, with linear fits overlaid. Downward black arrows 
mark the time of ejection of each component with the respective 
error bar. The gray vertical stripe indicates the epoch of the 7 -ray 
flare. 


2011 2012 2013 2014 

11' ' ' ' I ' ' 1' I ' ' ' ' II ' jj ' I ' ' ' ' I ' 

iN1 TN2 tN3 ^N4 

CO SI 
N2 
N1 N3 



Fig. 8 .— Light curves of 43 GHz model-fit components. Down¬ 
ward arrows and the gray vertical stripe indicate the same as in 
Figure 

properties are similar to those expected and observed 
previously for trailing components (Agudo et al. 2001; 
Gomez et al. 2001; Jorstad et al. 2005). Linear fits of 
separation versus time have been obtained for the other 
moving components, Nl, N2, N3, and N4, based on only 
those initial epochs at which an accurate position is ob¬ 
tained (see also Figure [^. This yields the estimates for 
the apparent velocities and times of ejection (epoch at 
which the component coincides with the core) listed in 

Tables 

Since we cannot directly measure the radial velocities 
of the jet features, a common approach to disentangle 
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the contributions of the component’s Lorentz factor and 
viewing angle in the observed proper motion is the use 
of the flux variability (e.g., Jorstad et al. 2005; Hov- 
atta et al. 2009). Following Jorstad et al. (2005), we 
use causality arguments to infer the variability Doppler 
factor 


_ sDl 
C^tvar(l “b '^) 


( 3 ) 


where s is the disk-equivalent angular diameter (where 
s=1.6o for a Gaussian component fit with FWHM=a 
measured at the epoch of maximum flux), and is 
the luminosity distance. The variability timescale is de¬ 
fined as ZXtvar—dt/ln(5'niax/*5’min)5 where Sniax and iS’niin 
are the measured maximum and minimum flux densi¬ 
ties, respectively, and dt is the time in years between 
'S'max and S'lnin (Burbidge et al. 1974). This defini¬ 


tion of Jvar 
density variability timescale corresponds 


is valid under the assumption that the flux 

to the light- 


travel time across the component, which is valid as long 
as the radiative cooling time is shorter than the light 
crossing time and expansion time. Combining the esti¬ 
mated value of Jvar with the measured apparent velocity, 
Papp = sin 0/(1 — f3 cos 6), where 9 and /3 are the view¬ 
ing angle and velocity (in units of the speed of light) of 
the component, we can calculate the variability Lorentz 
factor, Fvarj and viewing angle, 0var, using (Hovatta et 
al. 2009) 


r — 

-L va r — 




2 

app 


( 5 ^ 

^var 


2Jv 


and 


0var = arctan ■ 


2 / 3 , 


app 


app 


_ 1 ■ 


( 4 ) 


( 5 ) 


Physical parameters of the moving components ob¬ 
tained from this method are reported in Table 


TABLE 4 

Physical Parameters of Moving Jet Features 


Name 

^^var 

(yr) 

(masj 

^var 

^var 

n 

r var 

N1 

0.70 

0.14 

14.6 

3.9 

14.9 

N2 

1.12 

0.33 

22.4 

2.5 

19.6 

N3 

0.28 

0.09 

26.1 

2.2 

26.2 

N4 

0.20 

0.08 

30.3 

1.2 

17.3 


‘‘FWHM of the model-fit component calculated at the epoch of 
maximum flux 


4.2. Kinematics and Flux Density Variability 

By inspecting the light curves in Figurej^we can iden¬ 
tify two flaring periods in the core: a prolonged first 
event that extends from mid-2007 to the beginning of 
2009, and a second one between mid-2012 and the begin¬ 
ning of 2013, in coincidence with the main y-ray flare. 

The peak flux of the first flare occurs between 2008 
June and July, when both components CO and Cl in¬ 
crease their flux densities, reaching a combined value of 
^ 4.2 Jy. Due to the proximity of Cl to the core, it 
is not always possible for the model fitting routine to 
clearly distinguish the two components, leading to high 
uncertainties in the flux ratio of the two features, as well 


as uncertainties in the position of Cl. Because of this. 
Figure]^ also shows the combined flux density of the core 
and component Cl, providing the data needed to follow 
the total flux density within the core region of CTA 102. 

The second flare in the mm-wave core begins in mid- 
2012, reaching its peak flux density at the end of 2012 
October, close to the y-ray flare (see Figure ^ section 
[^. After the peak, the core region (CO plus Cl) remains 
m a high flux state until the last observing epoch, with 
a combined flux density oscillating around ^ 2 Jy. 

Both flares in the core region are associated with the 
appearance of subsequent superluminal components. In 
the case of the first such flare, component N1 appears 
as a bright and well-defined feature that moves along 
the jet at 14.9 ± 0.2 c (see Figs. 7]|8 and Table |^. We 
also note that the ejection of component N1 corresponds 
to a change in the innermost structure of the compact 
jet, after which component Cl is no longer detected for 
almost two years. The second core flare leads to the 
ejection of component N4, which is significantly weaker 
than component N1 and has the slowest proper motion 
of the analyzed components (see FigslTpland Table [^. 

The values of the variability Doppler factors listea in 
Table correspond to a progressive increase with time, 
from 14.6 for component N1 to 30.3 for component N4. 
Previous estimations of the va riability Doppler fact or in 
CTA 102 r ange between 15.6 (iH ovatta et al. 2009) and 
22.3 ± 4.5 (Jorstad et al. 2005|), making JN4 the super- 
luminal knot with the highest Doppler factor to date. 
According to our analysis this unusually large value is 
due to a progressive re-orientation in the direction of 
ejection of knots, from 0var = 3.9° for component N1 
to 0var = 1-2° for N4, which travels almost directly along 
the line of sight. This change in the jet orientation is 
readily apparent when analyzing each component’s po¬ 
sition angle shortly after the time of ejection, as well as 
their subsequent trajectories, as shown in Figure 

This smaller viewing angle of the jet with respect to 
the observer during the second radio flare, which appears 
to last until the end of our VLB A dataset (June 2014), 
is also in agreement with the significant differences ob¬ 
served between the ejections of components N1 and N4. 
While component N1 is clearly identified in the jet as a 
bright (> 2 Jy) component soon after its ejection, most 
of the increase in the total flux density during the sec¬ 
ond radio flare appears to be associated with the core 
region (CO-l-Cl), with component N4 representing only 
a small fraction of the flare. The smaller viewing angle 
of the jet also leads to a more difficult identification of 
component N4, which is not clearly discerned from the 
core until May 2014 (2014.33). Further support for the 
re-orientation of the jet toward the observed is also ob¬ 
tained from the analysis of the polarization, discussed in 
section |6l 


5. CROSS-CORRELATION ANALYSIS 

To quantify the relationship among the light curves at 
the different wave bands, we perform a discrete cross¬ 
correlation analysis. The z-transfor med discrete corre- 
lation function (ZDCF) described by Alexander (1997) 
has been designed for unevenly sampled light curves, as 
in our case. We use the publicly available zdcf-vl.2 and 
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Fig. 9.— Trajectories {left panel) and position angles as a function of distance from the core {right panel) of the moving components in 
the jet. 



Log (days) 


Fig. 10.— Z-transformed discrete correlation function between 
optical and 7 -ray data (upper panel) and UV and 7 -ray data (lower 
panel). In each panel we report the time lag corr esponding to the 
correlation peak with its respective l-tr error (see |Alexander|2013| 
for more details). 

plike-v4-0.f90 programwith a minimum number of 11 

http://www.weizmann.ac.il/weizsites/tal/research/software/ 


points inside each bin, as recommended for a meaningful 
statistical interpretation. We compute the DCF between 
each pair of light curves, including data from 100 days 
before to 100 days after the main 7 -ray outburst. Time 
sampling of the light curves ranges from one day for the 
7 -ray data to tens of days in the case of some other wave¬ 
bands (see Figure [^. 

Figure [T0| displays the ZDCF analysis for the optical- 
7 (upper panel) and UV -7 (lower panel) data. We find 
that the correlation peaks between the 7 -ray light curve 
and the optical and UV light curves give a time lag of 
0.70±1 and days, respectively, where a positive 

lag means that the 7 -ray variations lead. We therefore 
conclude that the variations at the three wavelengths are 
essentially coincident within the uncertainties. 

The sparser sampling of the X-ray data, as well as its 
double-peaked structure, precludes a reliable ZDCF anal¬ 
ysis. However, we note that the first X-ray data peak is 
coincident with the 7 -ray flare, and the second brighter 
X-ray flare occurs ^50 days later (see Figure]^. The 
triple-flare structure of the NIR light curve during the 
7 -ray flare (see Figure also prevents a unique inter¬ 
pretation of a cross-correlation analysis. Nevertheless, 
from inspection of the light curves, we see that the first, 
brightest peak in the NIR light curve (56193 MJD) is 
simultaneous with the 7 -ray outburst within an uncer¬ 
tainty of 1 day, corresponding to the time sampling of 
both light curves. 

We obtain no significant correlation between the 
millimeter-wave and 7 -ray light curves. This can be due 
to the different timescales associated with the emission at 
these wave bands, as suggested also for other blazars (e.g., 
1156-1-295, |Ramakrishnan et al.||2014|). The rise time for 
the millimeter wave band is of the order of months, while 
for the 7 -rays it is of the order of few days. We note, 
however, that the 1, 3, and 7 mm light curves contain a 
significant flare coincident with the 7 -ray outburst. 

6 . POLARIZED EMISSION 

Figure pi] shows the optical and millimeter-wave linear 
polarization between MJD 54000 and 56900, covering the 













































11 



Fig. 11. — Optical and millimeter-wave linear polarization over 
the period MJD 54000-56800. The first two panels display the de¬ 
gree of mm-wave and optical polarization, respectively. The third 
and fourth panels display the EVPAs at millimeter and optical 
wavelengths, respectively. The gray vertical stripe indicates the 
epoch of the main y-ray flare. 

period of the 7 -ray flare. To solve for the zbnTr ambiguity 
in the electric vector position angle (EVPA), we assume 
the slowest possible variation in time, applying a ±7r ro¬ 
tation between two consecutive measurements when the 
magnitude of the EVPA change would otherwise exceed 
7r/2. 

No significant increase in the degree of polarization at 
millimeter wavelengths is observed during the 7 -ray flare, 
but the EVPAs display a progressive rotation starting 
about one year prior to the 7 -ray flare. Figure 11 shows 
that between 2007 and mid-2011 the EVPAs at 3 and 
1 mm are distributed around a mean value of ~100°. Af¬ 
ter this, the polarization at millimeter wavelengths starts 
a slow rotation by almost 80° in one year (from 2011 
July to 2012 August), until the flare epoch. In coinci¬ 
dence with the 7 -ray flare, the rate of EVPA rotation in 
the VLBI core and stationary component Cl increases 
significantly, leading to a rotation of almost 200 ° in one 
year. Subsequently, component N4 appears and the EV¬ 
PAs of both Cl and N4 rotate again toward values similar 
to those at 1 and 3 mm, reaching ~ —100°. It is pos¬ 
sible that during the flare, while the new superluminal 
component N4 is crossing the core zone, the EVPAs of 
the innermost region at 7 mm rotate due to the passage 
of the component. After this, when N4 can be distin¬ 
guished from Cl and CO, the EVPAs at 7 mm again 
follow the general behavior of the EVPAs at shorter mil¬ 
limeter wavelengths. A similar discrepancy between the 
1-3 mm and 7 mm EVPAs occurs between mid-2009 and 
mid-2010, when component N1 is ejected and becomes 
brighter than the core until mid-2010 (see Figure!^. 

We can then distinguish the rapid rotation of polariza¬ 
tion vectors observed in the VLBI components at 7 mm 
after the flare from the slower rotation observed at 1 and 
3 mm. The latter leads to a rotation of the mm-EVPAs 



Fig. 12. — The upper panel displays, from top to bottom, the 
light-curve, degree and time evolution of polarization, and EVPA 
at optical frequencies. Each colored mark corresponds to the period 
over which we plot the U and Q Stokes parameters in the lower 
panel. Blue points mark a clockwise rotation cycle that occurs in 
coincidence with the total intensity peak. The gray vertical stripe 
indicates the time range of the main gamma-ray flare. 


of ^200° over 3 years (from 2011 July to 2014 August). 
This progressive rotation in the EVPAs can be produced 
by a change in the orientation of the innermost jet, which 
would be in agreement with the larger Doppler factor and 
smaller viewing angle of component N4 associated with 
the 7 -ray event, as discussed previously (see §|^. 

The optical polarization executes rapid and pro¬ 
nounced changes in both degree of polarizati on and 
EVPA associated with the 7 -ray flare. Figure dis¬ 
plays an expanded view of the optical polarization data 
near the time of the 7 -ray flare, with four different time 
ranges marked in different colors. Before the peak at op¬ 
tical frequencies on MJD 56194, the source undergoes a 
period of rapid changes in both total and polarized emis¬ 
sion (marked in red), and the EVPAs rotate by almost 
30°. The plot of Stokes parameters U vs. Q in Figure [T^ 
reveals a clear clockwise rotation of the EVPAs (marked 
in blue), in coincidence with the main flare in total flux 
and a rapid change in the degree of polarization. This 
clockwise rotation has been previously reported by |Lar- 


ionov et al. (2013a). 

It we assume a model in which a relativistic shock does 
not cover the entire cross-section of the jet and is moving 
down the jet following helical magnetic field lines, which 
also propagate downstream, then we expect to observe 
a rotation in the EVPA. This should be accompanied 
by a change in the degree of polarization, with a min¬ 
imum in the middle of the rotation, where the flaring 


/ty (e. 


Vlahakis 

2006 

Marscher et al. 

2008 Larionov 

et al. 

2 UJ 

L3b). Evidences ot a helical magnetic field in 


five circular polarization (Gabuzda et al.|[2008) and in a 
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gradient in the rotation me asure across the .jet width at 
about 7 mas from the core (Hovatta et al. | 2 dl 2 [ ), both in 
MOJAVE observations. 


7. SPECTRAL ENERGY DISTRIBUTION 

We have computed the spectral energy distribution 
(SED) of the source from millimeter to 7 -ray wavelengths 
at several epochs (see Figure 131: two epochs between 
MJD 54720 and MJD 55055, corresponding to 7 -ray qui¬ 
escent states; the epochs of the first 7 -ray flare (MJD 
55715-55721), the main flare (MJD 56193), and the third 
7 -ray flare (MJD 56392-56398); one epoch of a quiescent 
state between the first and second flares (MJD 56040- 
56047), and a second one between the second and third 
flares (MJD 56273-56280). For the main 7 -ray flare, 
all data are simultaneous except for the mm-wave data, 
which corresponds to MJD 56208. For the other epochs, 
we have considered a range of time (as indicated in Fig¬ 
ure 131 in order to cover the entire energy range. 



log V [Hz] 


Fig. 13.— Spectral energy distribution of CTA 102 from 
millimeter-wave to 7 -ray frequencies during the brightest 7 -ray 
flare (black) and at other observing epochs, as labeled in the figure 
and discussed in the text. 


By examining Figure |13[ we observe that, during the 
multi-wavelength flare in 2012 (black points) both, the 
synchrotron and the inverse Compton peaks increased. 
The synchrotron peak frequency during the flare is close 
to 10^^ Hz, at the upper end of the frequency range 
lO^^-lO^ '^ Hz of synchrotron peaks ob served in luminous 
blazars ( Ghisellini fc Tavecchio|2008 ). At the same time, 
we also observe a shift in the inverse Compton peak to 
higher frequencies, which leads to a hardening of the 
spectrum between 0.1 and 1 GeV. We note that neither 
of the other two, weaker 7 -ray flares display a similar 
shift in the peak of the inverse Compton spectrum. In 
particular, the weaker first and third 7 -ray flares dis¬ 
play a very similar inverse Compton spectrum, peaking 
at nearly the same frequency, with only a minor hard¬ 
ening toward higher energies during the third flare. The 
shift of the 7 -ray peak toward higher frequencies during 
the main 7 -ray flare can be explained by a change in the 
viewing angle - leading to an increase in Doppler factor 
- of the emitting region, in support of our hypothesis of 
a reorientation of the jet toward the line of sight during 
the multi-wavelength flare (see §§|^and|^. 


The ratio of the inverse Compton to synchrotron peaks 
is of the order of 10 , which is not sufficiently large to 
rule out synchrotron self-Compton scattering as the main 
mechanism for the production of the 7 -ray emission (e.g., 
Sikora et al.||2009 |. 


8 . DISCUSSION AND RESULTS 

We have presented a multi-wavelength polarimetric 
study of the quasar CTA 102 during an unprecedented 
7 -ray outburst that was observed between 2012 Septem¬ 
ber 23 and October 2. We find that the 7 -ray out¬ 
burst occurred simultaneously with flares from millime¬ 
ter to X-ray wavelengths, with the exception that the 
1 and 7 mm light curves peak almost one month af¬ 
ter the shorter-wavelength flares. However, all of the 
millimeter-wave light curves begin to increase before the 
shorter-wavelength outburst, but took longer times to 
reach maximum flux and then to decay. 

Our DCF analysis confirms the coincidence between 
the 7 -ray flare and the optical-UV flare. The same anal¬ 
ysis does not provide unambiguous correlation between 
the X-ray or NIR and the 7 -ray light curves because of 
the multi-peak structure of the flare at these two frequen¬ 
cies and the relatively sparse sampling. 

The Fermi LAT daily light curve reveals two more 
7 -ray outbursts apart from that in 2012 : one in 2011 
June and the other in 2013 April. Both outbursts are 
weaker and “orphan”. Only the bright outburst in 2012 
September-October is coincident with flares at the other 
wavelengths and with the emergence of a new superlu¬ 
minal knot from the radio core. 

We have combined our multi-wavelength study of light 
curves with an analysis of multi-epoch VLBA observa¬ 
tions at 43 GHz that provide the necessary angular res¬ 
olution to follow the evolution of the jet during the out¬ 
burst. In the 43 GHz VLBA images, we observe the 
ejection of multiple superluminal knots from the radio 
core during the analyzed period, but only one of these 
knots, N4, is associated with a 7 -ray flare. Component 
N4 was ejected in 2012.49 ± 0.11, within a time range 
between 47 and 127 days before the main 7 -ray flare in 
2012 (2012.73), when the radio core started to increase in 
flux density. The interaction between a traveling feature 
and the stationary radio core appears to ha ve triggered 


nan et al. 

2014 

Morozova et al. 2014 

1 and radiogalax- 

ies OGrandi et al. 

2012 Casadio et al. 

20151. However, 

not every ejection ot a new knot leads to a 7 -ray dare. 


For instance, it is not clear why a 7 -ray flare is associ¬ 
ated with N4 in CTA 102 and not with the other moving 
radio components. 

From the analysis of model-fit components at 43 GHz, 
we deduce that the jet changed its orientation with re¬ 
spect to the observer when component N4 was ejected. 
This is derived from an analysis of the variability Doppler 
factor and viewing angle, which indicates that a progres¬ 
sive increase in the Doppler factor occurred, caused by a 
re-orientation of the jet toward the line of sight. This led 
to a minimum viewing angle of 0 ~ 1 . 2 ° when component 
N4 was ejected during the 7 -ray outburst. This change 
in the orientation of the jet is supported by the observed 
progressive, slow rotation of the millimeter-wave EVPAs 
starting almost one year before the ejection of N4 and the 
7 -ray flare. We therefore conclude that the 7 -ray emis- 
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sion in CTA 102 is related to a decrease in the viewing 
angle of the jet. 

This correlation between 7 -ray activity and orientation 


iects 

Marscher et al. 

20081 Larionov et al. 20101 IRani 

et al. 

2L)14D, quasars 

Abdo et al. 

2010b Raiteri et al. 

mTf 

Jorstad et al. 

2013), and radio galaxies (ICasadio 

et al. 

2U15), although 

there are ditterent interpretations 


authors consider a bent or processing jet, while others 
suggest a helical jet with the radiating component follow¬ 
ing this helical path. A helical trajectory could also be 
the consequence of magn etic field lines t wisting around 
a conical or parabolic jet ( Vlahakis||2006 ). 

In the case of CTA 102, th ere are indications of a he¬ 
lical magri etic field structure (Gabuzda et al.|2008 Hov- 


atta et aL]|2012 1 . We associate the fast variability in the 


polarized optical emission, as well as the clockwise rota¬ 
tion displayed in the EVPAs during the outburst, with 
the helical path followed by the superluminal component 
in its motion along the outwardly propagating magnetic 
field lines. On the other hand, a number of similar rota¬ 
tions of the mm-wave and optical polarization vectors oc¬ 
curred in both the clockwise and counterclockwise direc¬ 
tions over the entire 2004-2014 monitoring period. This 


lent magnetic field ( 

Jones|1988 

D’Arcangelo et al. 

2007 

Marscher 

2014 

1. Early results trom the KoboPo 

pro- 


7 -ray flares c an be produced by a random walk process, 
some are not (Blinov et al.||2015). If the rotation associ¬ 
ated with a 7 -ray flare is caused by a helical geometry of 
the magnetic field, then future such outbursts should be 
accompanied by similar clockwise rotations. 

The observed long-term rotation in millimeter-wave 
polarization vectors, together with the slower proper mo¬ 
tion associated with component N4, suggest a change 
in the jet orientation, so that it becomes more closely 
aligned with the line of sight during the ejection of com¬ 
ponent N4 and the multi-wavelength flare. 

The close timing of the 7 -ray, X-ray, UV, and optical 
flares suggests co-spatiality of the emission at all these 
frequencies. Knot N4 was 0.025 to 0.07 mas downstream 
of the core when the 7 -ray flare occurred, i.e., it had not 
yet reached feature Cl at ^0.1 mas. This is confirmed 
by the increase in flux density in the 7 mm core during 
the 7 -ray outburst. Hence we conclude that the bright 
7 -ray outburst occurred inside the mm-wave core region. 

We observe component N4 for the first time in the 
VLB A images on 2013 April (MJD 56398), when it was 
located at r ~0.12 mas. The 7 -ray flare in 2013 April oc¬ 
curred between MJD 56387 and 56394. Therefore, a pos¬ 
sible interpretation of this flare is t he passage of compo- 
nent N4 through Cl, interpreted by Fromm et al. (2013aI 
as a possible recollimation shock. 

If the radio core were located within ^ 1 pc of the 
black hole (BH), then the accretion disk or the broad 
line region could provide the necessary photon field to 
explain the high energy emission through external Comp¬ 
ton scattering. The 43 GHz radio core in CTA 102 must 
be coincident with, or downstream of, the 86 GHz core 
that is located at a distance of 7.5 ± 3.2 pc (~ 8.5 x 10"^ 
gravitational radii for a BH mass of ~ 8.5 x IO^Mq; fZa-| 


maninasab et al. 20141 from the BH ( Fromm et al.|20i5 ). 
A similar scaled distance, ~ 10"*^-10® gravitational radii, 
has been de termined also for two radio galaxies, 3C 111 
and 3 C 120 (Marscher et al.|2002 Chatterjee et al.|2009 


201^ and two blazars, BL Lac and 3C 279 (Marscher 


et al. 


gle o: 


20081 lAbdo et al.||2010b|) . For a mean viewing an- 


the jet of C'lA 102 of 2.6° (Jorstad et al. 2005 


Fromm et al. |2015[ ), the distance of 1N4 from CO is 4.6-1. 


pc, hence the 7 -ray outburst took place more than 12 pc 
from the BH. At this location, there should be a negli¬ 
gible contribution of photons from the disk or the broad 
line region, nor from the dusty t orus (located ^ 1.6 pc 
from the BH; Pacciani et al.|2014), for external Compton 
scattering to produce the high energy flare. The lack of 
a suitably strong external source of photons favors syn¬ 
chrotron self-Compton (SSC) scattering of NIR to UV 
photons by electrons in the jet with energies ~ 10 times 
the rest-mass energy, as the source of the 7 -ray emission. 
The ratio of 7 -ray to infrared (synchrotron) luminosity 
is < 10, sufficiently low to be consistent with the SSC 
process. 


9. CONCLUSION 

Our study of the time variability of the multi¬ 
wavelength flux and linear polarization of the quasar 
CTA 102 confirms its erratic blazar nature, revealing 
both strong connections across wave bands in one out¬ 
burst and no obvious connections for other events. The 
bright 7 -ray outburst in late 2012 was accompanied by 
contemporaneous flares at longer wavelengths up to at 
least 8 mm, with the increase in mm-wave flux starting 
before the 7 -ray activity. The polarization vector at both 
optical and millimeter wavelengths rotated from the time 
of the 7 -ray peak until ~ 150 days later. A new superlu- 
minally moving knot, N4 — the feature with the highest 
Doppler beaming factor during our monitoring, accord¬ 
ing to our analysis — was coincident with the core in the 
43 GHz VLBA images 47-127 days prior to the 7 -ray 
peak. We conclude that the outburst was so luminous 
because the jet (or, at least, the portion of the jet where 
most of the emission occurs) had shifted to a direction 
closer to the line of sight than was previously the case. 
The time delay between the epoch when N4 crossed the 
centroid of the core (feature GO) and the epoch of peak 7 - 
ray emission implies that the main flare took place > 12 
pc from the black hole. At this distance, the only plausi¬ 
ble source of seed photons for inverse Gompton scattering 
is NIR to UV emission from the jet itself. The ratio of 
7 -ray to infrared luminosity is only ~ 10 at the peak 
of the outburst, low enough to be consistent with SSG 
high-energy emission. 

Multiple superluminal knots appeared in the jet dur¬ 
ing the 7 years covered by our VLBA observations. 
These include a very bright component (Nl) ejected in 
2009.12±0.02 and associated with a significant mm-wave 
flare in the core region (~4.2 Jy). Yet only component 
N4 is related to a flare at 7 -ray energies. Two strong “or¬ 
phan” 7 -ray flares have no apparent optical counterparts. 
A strong mm-wave event with neither a 7 -ray nor optical 
counterpart can be explained by an inability of the event 
to accelerate electrons up to energies ~ 10 ^ mc^ needed 
to radiate at such frequencies, although the reason for 
this inability is unclear. Orphan 7 -ray flares might be 
explained by a knot crossing a region where there is a 
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higher loca l density of seed photons for inverse Compton 
scattering ( Marscher et al.|2010 MacDonald et al.|2015 ). 
Indeed, the second orphan hare corresponds to the tirrie 
of passage of knot N4 through stationary feature Cl (lo¬ 
cated ^0.1 mas from the core), which could be such a 
region. 

During the multi-wavelength outburst we observe 
intra-day variability in the optical polarized emission, as 
well as a clockwise rotation in optical EVPAs. This ro¬ 
tation could be caused by a spiral path traced by the 
knot moving along helical magnetic field lines that prop¬ 
agate outwards relativistically. Alternatively, the various 
rotations of the polarization vector seen in our dataset, 
which are in the clockwise and counterclockwise direc¬ 
tion over different time ranges, could be mainly random 
walks caused by a turbulent magnetic field. 

CTA 102 displays the complex behavior characteristic 
of the blazar class of active galactic nuclei. Nevertheless, 
we have found possible connections between variations 
in the multi-wavelength flux and polarization and in the 
structure of the jet in some events. Continued monitoring 
of CTA 102 and other bright blazars at multiple wave 
bands with as dense sampling as possible, combined with 
mm-wave VLBI imaging, can eventually determine which 
connections are robust and the extent to which stochastic 
processes dominate the behavior of blazars. 
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